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Abstract 

Expansions of the gravitational field arising from the development of asymptotically Eu¬ 
clidean, time symmetric, conformally flat initial data are calculated in a neighbourhood of 
spatial and null infinities up to order 6. To this ends a certain representation of spatial 
infinity as a cylinder is used. This set up is based on the properties of conformal geodesics. 
It is found that these expansions suggest that null infinity has to be non-smooth unless the 
Newman-Penrose constants of the spacetime, and some other higher order quantities of the 
spacetime vanish. As a consequence of these results it is conjectured that similar conditions 
occur if one were to take the expansions to even higher orders. Furthermore, the smoothness 
conditions obtained suggest that if a time symmetric initial data which is conformally flat in 
a neighbourhood of spatial infinity yields a smooth null infinity, then the initial data must 
in fact be Schwarzschildean around spatial infinity. 


1 Introduction 


Penrose introduced the seminal idea that the gravitational field of isolated systems can be conve¬ 
niently described by means of the notion of asymptotic simplicity HU Q Central to the concept 
of asymptotic simplicity is the idea —expectation— that the conformal boundary of the space- 
time null infinity, — should possess a smooth differentiable structure. This approach to the 
description of isolated bodies in General Relativity is usually known as Penrose’s proposal —see 
e.g. 118, 19 . Static and stationary spacetimes have been shown to be (weakly) asymptotically 


simple, with a smooth null infinity jj]. However, the main purpose behind introducing the con¬ 
cept of asymptotic simplicity is to provide a suitable framework for the discussion of radiation. 
In spite of its elegance and aesthetical appeal, Penrose’s proposal is of little use if one is not 
able to prove that there exists a big family of non-trivial (in the radiative sense) asymptotically 
simple solutions to the Einstein Field Equations. A programme aimed to investigate the existence 
of such solutions, and provide a conclusive “answer” to Penrose’s proposal has been started by 
Friedrich —see e.g. [jn], |l2], [lil bl [bil - [T?I . His strategy is based on the use of the so-called 
Conformal Field Equations, which allow us to work and prove existence statements directly in 
the conformally rescaled, “unphysical” spacetime. Along this lines, Friedrich has been able to 
prove a semiglobal existence result that ensures that hyperboloidal initial data close to Minkowski 
data yields an asymptotically simple development which includes the point i + (future timelike 
infinity) jn|. Recently, Chrusciel & Delay Q] -retaking an old idea by Cutler & Wald Q 
using a refined version of some initial data set constructed by Corvino have been able to prove 
the existence of a big class of non-trivial (radiative) asymptotically simple spacetimes. Corvino’s 

* E-mail address: jav@aei-potsdam.mpg.de 

1 It will be assumed that the reader is familiar with the ideas of the so-called conformal framework to describe 
the properties of isolated bodies and the concept of asymptotic flatness. For a recent review, the reader is remitted 
to 
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initial data is constructed so that it is Schwarzschildean in a neighbourhood of spatial infinity, i°. 
This means that the radiation content in the spacetime arising from the development of the data 
is somehow special. This can be seen directly from the fact that the Newman-Penrose constants 
of the spacetime are zero Bizi- Clrrusciel & Delay’s result is no doubt very important. However, 
it is not as general as one would like. 

It has been suspected for a long time now that the region of spacetime where null infinity 
and spatial infinity meet is somehow problematic —see e.g. ^3|. From the analysis of the 
hyperboloidal initial value problem it turns out that the smoothness of null infinity is preserved 
by the evolution if smooth data sufficiently close to Minkowski data are prescribed. The latter 
indicates that somehow the “decision” of having a smooth structure at null infinity is made in 
an arbitrarily small neighbourhood of spatial infinity. In some sense, Corvino’s data avoids all 
the intricacies and complications of this region of spacetime by setting the asymptotic end in the 
simplest way which is consistent with the presence of a non-vanishing ADM mass on the initial 
hypersurface. In connection with this, Friedrich |jl^| has performed a detailed first analysis of 
the behaviour of the gravitational field arising from asymptotically Euclidean, time symmetric 
initial data in the region where null infinity “touches” spatial infinity. By means of a novel 
representation of spatial infinity in which the point i° of the standard conformal picture is blown 
up to a cylinder I — the cylinder at spatial infinity— a certain regularity condition on the 
Cotton-Bach tensor and symmetrised higher order derivatives of it has been obtained. The hope 
was that this regularity condition would ensure the smoothness of null infinity, at least in the 
region close to i° and J'. A subsequent analysis by Friedrich & Kannar 1201 of the first orders of 
some expansions that can be obtained by evaluating the Conformal Field Equations at I, lead to 
conjecture that Friedrich’s regularity condition is the only condition one has to impose on time 
symmetric, initial data possessing an analytic compactification in order to obtain a development 
with smooth null infinity JFIl . More precisely, 


Conjecture (Friedrich, 2002). There exists an integer A;* > 0 such hat for given k < fc* 
the time evolution of an asymptotically Euclidean, time symmetric, conformally smooth initial 
data set admits a conformal extension to null infinity of class C k near spacelike infinity, if the 
Cotton-Bach .spinorf] satisfies the condition, 


D(a s b s ' ' ' D a - l b 1 b a bcd)(i') S 0, 1, . . . 

for a ceHain integer s* = s*(fc). If the extension is of class C°° then the condition should hold 
to all orders. 


The objective of this paper is to provide a further insight into the conjecture above. It will turn 
out that the conjecture, as it stands, is false. In order to see why this is the case, the expansions of 
Kannar & Friedrich will be carried to an even higher order. This requires the implementation of 
the Conformal Field Equations on a computer algebra system (Maple V). It should be emphasized 
that despite the use of the computer to perform the expansions, the results here presented are 
exact up to the order carried out. In order to simplify our discussion, the analysis will be restricted 
to developments of time symmetric conformally flat initial data. The time symmetry requirement 
stems from the fact that Friedrich’s analysis has only yet been carried for this class of data. A 
similar analysis of initial data with non-vanishing second fundamental form lies still in the future. 
However, some first steps have already been carried out j8|. The conformal flatness of the data 
ensures that the initial data satisfies the regularity condition trivially. Again, the construction 
of non-conformally flat data satisfying Friedrich’s regularity condition is a non-trivial endeavour 
whose undertaking will be left for future studies. In the light of the results here presented, it turns 
out that the restriction to the class of conformally flat data is not a drawback. Furthermore, it 
is not hard to guess how the results could generalise in the case of general time symmetric data. 

The principal result of our investigation is the following, 

2 The Cotton-Bach tensor Bijk is related to the Cotton-Bach spinor via: 

t^ijk 1 ' babce^df T t> a bdf^ce- 

This correspondence is carried out by the Infeld-van der Waerden symbols. 
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Theorem (Main theorem). Necessary conditions for the development of initial data which are 
time symmetric, conformally flat in a neighbourhood B a {i ) of (spatial) infinity to be smooth at 
the intersection of null infinity and spatial infinity are that the Newman-Penrose constants G k , 
k = 0,... ,4 and the higher order Newman-Penrose constants, G^, k = 0,... ,6 vanish. 

A more precise formulation of the theorem, including the definition of the Newman-Penrose 
constants and the higher order Newman-Penrose constants in terms of the initial data will be 
given in the main text. 

We note that the previous conditions are fulfilled by the Schwarzschild initial data. The 
theorem has been obtained from the analysis, up to order p = 6 of expansions constructed 
from the solutions of the transport equations induced by the conformal field equations upon 
evaluation on the cylinder at spatial infinity. From the evidence provided by the equations it is 
not unreasonable to conjecture that similar conditions arise if one were to obtain expansions to 
even higher orders. In 0 it was shown that the whole set up of the cylinder at spatial infinity 
is completely regular for Schwarzschildean data. Thus, the hypothetic new conditions must be 
satisfied by the Schwarzschild initial data. It is on this ground that the following conjecture is 
put forward: 

Conjecture (New conjecture). If an initial data set which is time symmetric and conformally 
flat in a neighbourhood B a (i) of the point i yields a development with a smooth null infinity, then 
the initial data are in fact Schwarzschildean in B a {i). 

Again, a more technical version of the conjecture is given in the main text. 

The article is structured as follows: in section 2 a brief summary of the description of spacetime 
in the neighbourhood of spatial and null infinities in terms of the cylinder at spatial infinity is 
given. This digest has the intention of providing all the tools needed for the calculations described 
in the present article. Particular attention is paid to the spatial 2-spinor formalism and to the 
expansions of functions on S 3 in terms of unitary representations of SU{ 2, C). The reader is, 
in any case, remitted to jl^] for a more extensive discussion. In section 3, the conformal field 
equations written in the conformal geodesic gauge are discussed. The initial data for the latter 
in the case of an asymptotically Euclidean, time symmetric, conformally flat initial hypersurface 
are described. The transport equations implied by the conformal field equations on the cylinder 
at spatial infinity are introduced in section 4. In order to properly motivate the calculations 
carried out in later sections, a careful discussion of their properties is undertaken. Section 5 
contains the new results to be presented in the article. Here a description of the solution of the 
transport equations is done. Due to the large size of the expressions involved, the description will 
be focused on what we believe are the most relevant features of the solutions. However, it should 
be emphasized that everything has been explicitly calculated. The main conclusions extracted 
from the calculations is presented as our main theorem. In order to understand the meaning 
of the conditions presented in the main theorem, the Schwarzschild solution is discussed in this 
context, and a conjecture is formulated. Finally, in section 6, some conclusions and extensions of 
the present work are considered. There is also an appendix in which the implementation of the 
transport equations and how to solve them in the computer algebra system Maple V is briefly 
considered. 


2 Spacetime in a neighbourhood of spatial infinity 

Let (M, g kV ) be a vacuum spacetime arising as the development of asymptotically Euclidean, time 
symmetric, conformally analytic initial data (S,h at p). Later, we will restrict the class of initial 
data sets under discussion to those conformally flat around infinity. Assume for simplicity that S 
possesses only one asymptotic end. Let i be the infinity corresponding to that end. The point i 
is obtained by conformally compactifying the initial hypersurface S , with an analytic conformal 
factor f l which can be obtained from solving the time symmetric constraint equations. The 
compact 3 dimensional manifold obtained in this way will be denoted by S , and the conformally 
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rescaled 3-metric by h a p. Assume that the 3-metric h a p is analytic in an open ball B a (i) of radius 
a centered on i. Let p denote the geodesic distance along geodesics starting at i. The radius a of 
the ball B a (i) is chosen such that B a (i) is geodesically convex. 

Now, let N D M be the domain of influence of the ball B a {i ) D S. Intuitively, one expects N 
to cover a region of spacetime “close to null and spatial infinities”. In reference ]l7j it has been 
shown that once the time symmetric constraint equations have been solved, a certain gauge based 
on the properties of conformal geodesics can be introduced. Let r be the parameter along these 
curves. This gauge has the property of producing a conformal factor 0 which can be in turn used 
to rescale the region N to obtain a “finite representation” N, of spacetime in a neighbourhood of 
spatial and null infinities. The relevant conformal factor is given by, 

e = « -ln (l -r2 ^) ’ W 

where ui is given by, 

u = — =, ( 2 ) 

and k is a smooth function such that k = n'p with k'(i) = 1. It contains the remaining piece of 
conformal freedom in our setting. 

Throughout this work, space 2-spinors will be systematically used. In order to avoid problems 
with vanishing frame vectors on surfaces diffeomorphic to spheres, our discussion will be carried 
not on N but on a subbundle C a ^ K of the frame bundle over N. The subbundle C a ^ K can be 
shown to be a 5-dimensional submanifold of 1 x 1 x SU( 2,C) with structure group 17(1). More 
precisely, we define C a , K to be given by, 


C a = I (t, p, t) £ R x R x 517 (2, C)| 0 < p < a, (3) 

l K K, ) 

The projection 7r' of C a , K into N corresponds to the Hopf map 517(2, C) —> 517(2, C)/U(l) ~ S 2 . 
Scalar fields and tensorial fields on N are lifted to C a ,K- Their “angular” dependence will be then 
given in terms of functions of t £ 517(2, C). 

The manifold C a ,K has the following important submanifolds, 

I = {(r,P,t)€C a , K \p = 0, | t | < 1 }, (4a) 

/=t = {(r,p,t) £ (7 aj re | p = 0, r = ±l}, (4b) 

= |(r,p,t) £ C Ql re | p > 0, r = ±^j. (4c) 


2.1 Space 2-spinors 

Consider the antisymmetric spinors e a b, £ ab , a,b = 0,1. These satisfy eoi = 1) e 01 = 1. Let r aa 
denote the tangent vector to the conformal geodesics parametrised by r. We set, 


c o c 0' 


(5) 


In order to define introduce differential operators on 517(2, C) we consider the basis 



( 6 ) 


of the Lie algebra of SU(2,C). Note that in particular, 113 generates 17(1). Denote by Zi, i = 1, 2, 3 
the real left invariant vector field generated by rq on SU(2,C). We define the, following complex 
vector fields 

X + = -(Z 2 + iZ 1 ), X_=-{Z 2 -iZ 1 ), X = -2 iZ 3 . (7) 

These are such that [X + ,A'_] = —X. With the help of X± one can construct the following 
(frame) vectors fields, 

Cab = C° aa , dr + cl a , dp + c+, A+ + c~ a , X-, (8) 
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on C a ,K- The use of a space spinor formalism based on the vector field \f2d T = T aa c aa ' allows 
to perform our whole discussion in terms of quantities without primed indices. Accordingly, we 
write 


Caa' 


^j^aa'dr T a ,C ab 


(9) 


with c ab — T ( a & c b )b' = c° ab d T +c 1 ab dp+c+ b X + +c ab X_. The connection is represented by coefficients 
T abed which can be decomposed in the form, 


T'abcd (jiabed X(ab)cd ) fc/od, 


( 10 ) 


where the fields entering in the decomposition possess the following symmetries: Xabcd = Xab{cd)i 
£,abcd = £(ab)(cd)> fab = f( a b)- The curvature will be described by the rescaled conformal Weyl 
spinor (j> a bcd = 4>(abcd ) j and by the spinor field Q a bcd = O a fc(cd) which encodes information relative 
to the Ricci part of the curvature. For the purpose of writing the field equations it will be 
customary to consider it decomposed in terms of Q g 9 cd and Q( ab ) cd . 

For latter use it is noted that an arbitrary four indices spinor X abcd can be written in terms of 
the “elementary spinors” e l abcd , e ac x bd + e bd x ac , e ac y bd + e bd y ac , e ac z bd + e bd z ac , and h abcd where, 

X a b = '^ e (a e b ) i Vab = ~y^ 6 a £ b > Z ab = ~^ £ a € b > (H) 

and, 

e abcd = e (o e b e c e d) > h abcd = — £ a (c e d)b- (12) 

The notation ( abcd)i means that the indices are to be symmetrised and then i of them set to 1. 
We write, 


Xabcd — X 0 e° abcd + X\ Cabcd + X 2 e 2 abcd + X 3 el bcd + X 4 e^ bcd 

T X x (c ac x bd T c bd x ac ) T Xy{c ac y bd -t- e b ^|/ ac ) T- X z {t ac z bd T c bd z ac ) 

+ X b h abcd . (13) 


2.2 Expansions of functions on C a , K 

In order to obtain our expansions of the gravitational field around spatial infinity, we will require 
to decompose the functions arising into their diverse spherical (harmonic) sectors. Any function 
real analytic complex value function / on 5/7(2, C) can be expanded in terms of some functions 
\/m + 1 T m k - forming a complete set in T 2 (^z, SU(2, C)) where /z is the standard Haar measure 
in 5/7(2,C). One has, 

oo m m 

f(t) = X] V 

m—0 j =0 k —0 

The functions T m k j can be shown to be related with the standard spin-weighted spherical harmon- 
ics. Using the fact that the group 5/7(2, C) is diffeonrorphic to 5 3 , one can use the coordinates 
(C, C > a ) t° write a given t € 5/7(2, C) as, 


1 

v/1 + CC 


e ia _ ie~ ia C \ 
ie ia C e ia ) ‘ 


(14) 


The latter yields the following correspondence rule between the functions T m 7 and the spin- 
weighted spherical harmonics s Yi n \ 


Y 

s 1 n/. 


H) 


s-\-2n—m 


2n T 1 
47T 


To. 


(15) 
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The operators X± and A' (related to the NP 9 and 9 operators) introduced in the previous 
sections can be seen to yield, upon application to the T m functions the following, 

X+T m k j =!3 mtj T m k j _ i, (16a) 

X - T m k j = -Pm,j+iT m k j+1 , with f3 m j = y/j(m- j + 1). (16b) 

A function / on SU(2 , C) is said to have spin weight s if X f = 2s f. This definition can be readily 
extended to functions on C ajK . As it will be seen later, all the quantities we will work with will 
have a well defined spin weight. Let / be an analytic function / on C a , K with an integer spin 
weight s. 

Now, consider a spinorial symmetric analytic function on C a , K , A 0i ... 02t . with essential com¬ 
ponents A j = A ( ai ,---a 2 r)j> 0 < j < 2r, of spin weight s = r — j. Then, the components of the 
function will possess expansions of the form, 


OO 

A j = 

p—0 


(17) 


where the coefficients A j }P can in turn be decomposed in terms of the functions as 

q(p) 2q 

-\j',P = ^i,p;2g,fc^2n q-r+j' (1®) 

q=\r-j\ fc=0 

where 0 < |?— j\ < g(p) < oo. An expansion of the latter type will be referred to as an expansion 
of type q{p). 

The conformal field equations are nonlinear. Thus, when expanding them, one finds products 
of T-functions. These products can in turn be reexpressed as a linear combination of T’s. More 
precisely, one has: 

Lemma 1. Multiplying T functions. The following holds, 

T 1 v T k2 _ 

2ni l\ 2 n 2 I 2 

ni+n 2 

(—1 ) n+ni+n 2 C(ni,m — h;n 2 , «2 — h m , n, m + n 2 — li — I 2 ) 

n=q 0 

x C(m,m — ki;n 2 ,n 2 — fe; n,m + n 2 — fci — fe) 

rp n+k 1 +k 2 -ni-n 2 /iq\ 

^ 2 n 71+Z1+Z2—ni —ri2 ’ 

where qo = max{\ni — U 2 \, ni + U 2 — ki — k 2 , ni + U 2 — h — 12 }, and C{l\,m\ \ I 2 , m 2 ; Z, m) denotes 
the standard Clebsch-Gordan coefficients of SU( 2,C) [] 


3 The conformal evolution equations 

Using the conformal geodesic gauge and the two spinor decomposition, it can be shown that the 
extended conformal field equations given in 111 imply the following evolution equations for the 

3 Some other used notations in the physics literature are: 

C(h,mi;l2,m2;l, m) =< li, h;mi,m,2\l, m >, 

= C(li,l2,l\mi,m2,m). 
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unknowns v = (d£ b , £, a bcd, fab, X{ab)cd , ©(ah)cd, 0/ cd ), A* = 0, 

drC° h = -X {ab) ef C° ef - fab, 
drCab = —X(, ^ 


= 0 , 1 .±. 


3b) C e/ 5 

e f 


1 


d T fabcd — X( a &) fefcd "T ^^{iacX(bd)ef T ^bdX(ac)ef)f 


ef 


^^X( a (,)( c fd)e y( e acQ/ b d + e bd®/ ac ) i&Pabcd, 


drfab ~ X( a b) ^ fef + / a b’ 

3rX(ab)cd = ~ X( a {,) Xefcd 0(ed)ab "t~ © Pabcd, 

&T®(ab)cd X( a b) ©( a ^) e / ^r^Vabcd T i's/^d ( a ftb)cde, 

® T ®g 9 ab = — X( a 6) ^® g 9 ef + * Wabef , 


.P) f 


where 


Vabcd — iff abed T "hi ^~b ^ 


<2 Wb'c'd')) Habcd — ^{ffabed T a T b TT - 


(20a) 

(20b) 

(20c) 

(20d) 

(20e) 

(20f) 

(20g) 


f’a'b'c'd'), (21) 


denoting respectively the electric and magnetic part of of 4> a bcd- The quantities 0, <9 T 0 and d a b, 
given by formulae (Ij) and ( ph| ) are known directly from the initial data. Thus, the equations 
(20a)-(20g) are essentially ordinary differential equations for the components of the vector v. 

The most important part of the propagation equations corresponds to the evolution equations 
for the spinor (f> a bcd derived from the Bianchi identities Bianchi propagation equations: 

(V2 - 2c° m )d T (j)o + 2cq 0 <9 t <(>i - 2co 1 <9 Q </> 0 + 2co O <9 Q 0i 
= ( 2 Toon — 8 Tioio)^o + (4Toooi + 8 Tiooo)<^i — 6 roooo 02 j (22a) 

Vzdrfii - Ciidrfj )0 + CQ 0 d T (p 2 ~ c“ 48^0 + CQ Q d a (p 2 

= — (4Tino + fn)4>o + (2Toon + 4Tnoo — 2/oi)0i + 3/oo</>2 — 2roooo03> (22b) 

V 2 d T (j )2 - c° n d T (l )i + c° 00 d T (t>z - i + Co O d a 0 3 

= —Tiin^o ^ 2 (Thoi + fii)<pi + 3(Tooii + Tnoo)02 — 2(Toooi — /oo)03 — Toooo04(,22c) 
V2dr4>3 -Ci 1 d T (j>2 + CQ 0 d T (p4 - Clidalj)2 + CQ Q d a (p4 
= — 2rnn</>i — 3fii(f>2 + (2Tnoo + 4Toon + 2 /oi) 0 3 — (4Toooi — /oo)04, (22d) 

(V2 + 2 cQ 1 )d T 4>4 - 2c° 11 d T (j) 3 + 2cq 1 9 q </> 4 - 2c“ 1 d Q 0 3 
= —6riin^>2 + (4Tino + 8r 0 m)(/>3 + (2Tnoo — 8roioi)04- (22e) 

To the latter we add a set of three equations, also implied by the Bianchi identities which we 
refer as to the Bianchi constraint equations, 

C n <9 T 0O 2 CqiC1t 01 T CQQtlr02 T c^d a (j)Q 2c m d a (f>i T c oo d a 0 2 
= —(2T ( 01 ) 11 — 4r mo )0o + (2Toon — 4T( 01)01 — 4Thoo)0i + 6 r( 0 i)oo <( > 2 — 2r O ooo0ij?3a) 

C?i<9 T 0i - 2 Co 1 <9 t </>2 + C° m d r (t> 3 + Ci^a^i - 2 Coi<9q,</>2 + c 30 d a (t) 3 
= Tnn^o — (4T( 0 i)ii — 2T lwl )cj) 1 + 3(r 0 on — Tnoo)02 — (2r 00 oi — 4r( 0 i)oo)</ > 3 
~T O ooo04> (23b) 

c?A 0 2 - 2 cQ 1 d T (/)3 + c° m d T <j)4 + c^dafa - 2 cg' 1 0 a 0 3 + c“ o S Q 0 4 
= 2T lin <j> 1 — 6T(oi)n02 + (4Toon + 4T( 0 i)oi — 2r 11O o)0 3 — (4r 0 ooi — 2T(oi)oo)04-(23c) 


3.1 The initial data 

As pointed out in the introduction, only asymptotically Euclidean, time symmetric, analytically 
conformally flat initial data will be considered in our discussion. A number of simplification 


7 




arise under these assumptions, 
hypersurface can be written as, 


In particular, around i, the conformal factor Q of the initial 


n = 


(1 + pW ) 2 ’ 


(24) 


where W{i ) = m/2, m the ADM mass of the initial hypersurface. The function W satisfies the 
Yamabe equation, which under our assumptions reduces to the Laplace equation. Therefore W 
is harmonic, and thus can be written as, 


oo 2 p 

" = y + ~\ w p' 2 p^ t, 2p pi? 

p= 1 fe=0 


(25) 


where the coefficients w P: 2 P ,k, P — 1, 2 ,..k = 0,..., 2p complex numbers satisfying the regularity 
condition r w p , 2 P ,k — (-i y +k w p . , 2 p, 2 p-k so that IT is a real valued function. 

As mentioned in section 2 , a crucial property of our set up based on the properties of conformal 
geodesics is that it renders a conformal factor 0 —see formula (il|) for the region of spacetime 
under discussion. Furthermore, solving the con form al ge odes ic equations also yields a 1-form 
dab , which appears in the propagation equations (201) and (|20^ ). Under our assumptions of time 
symmetry and conformal flatness it is given by, 


dab — 2 p 


%ab P B) ab 

(1 + pW) 3 


(26) 


Once t he fu nc tion n of section 2 has been chosen, the initial data for the conformal propagation 
equations ( | 20 a|) -( 20 g) is given by, 

J 2 


O abed — D^abl^cd)^! 


n 


tfiabed — ^2 D( a bD c d)Q'> 

Cab = 0 , c\ b = KX ab , 


C — Zabi 

P 


C — Vab i 

P 


K 1 

I 7> {^ac^bd T ^bd^ac) 71 (^ac-Dfcd/'C Y ^bdB ac tvj 
( Ip 2 k 


f,abcd — V^2 \ - 
X(ab)cd = 0 ; fab — D ab K. 

where D a b , the spinorial covariant derivative of the initial hypersurface S is given by, 

Babpcd — •Eab^pPed Y ^ab-^+Pcd Y JJab^- — Pcd 7afe cPed 7af> dPec^ 

P P 

where the flat connection coefficients 7 a b c d are given by, 


(27a) 

(27b) 

(27c) 

(27d) 

(27e) 

(27f) 

(28) 


7 abed — n (tac^bd Y Cfod^ac); 
2 P 


(29) 


for a given differentiable spinorial function p ab . 


4 The transport equations 

The equations ( |20a| )-( |20g| ) can be concisely written in the form, 

d T v = Kv Y Q(v, v) + L(f>, (30) 


where K and Q are respectively a linear and a quadratic function with constant coefficients, 
whereas L is a linear function depending on the coordinates via 0, <9 r 0 and d a b . For the Bianchi 
propagation equations ( 22 a)~(| 22 e|) one can write, 


V2 Ed T cl> + = B(T abcd )0, 

8 


(31) 







where now, E denotes the (5 x 5) unit matrix, A ab c p b are (5 x 5) matrices, and B(T abcd ) is a 
linear (5 x 5)-matrix valued function of the connection coefficients T abcd . Similarly, the Bianchi 
constraint equations (23a)-(E3q) can be written as, 


F ab c£ b d^ = H(T abcd ), 


(32) 


where now F ab c p b denotes a (3 x 5) matrix, and H(T abcd ) denotes a (3 x 5) matrix valued function 
of the connection. 

In the sequel, given an unknown u we will write, = u\i. The objects 0, d T 0 and d ab 
from which L is constructed, vanish on /. Thus, L^ = L\i, and consequently the equations (|30|) 
and (31) decouple from each other. The system of equations for the v unknowns, equation (|30|), 
turns out to be an interior system upon evaluation on the cylinder at spatial infinity. Its initial 
data can be read from the restriction of the initial data ( [27a|) -( 27f) to I. It can be seen that 
this restriction, irrespectively of the choice of the function k coincides with the initial data of 
Minkowski spacetime. With this information in hand, the system for can be readily solved 
yielding, 


y(°) = o 

^•( ab)cd ’ 


f (°) — o — 0 

Jab ~ u ’ ^abcd ~ u ’ 


® { l d = 

( C ah) (0) = -TXab, (Cl b r = 0, (C“ h ) (0) = Vab, (c+ h ) (0) = Z ab . 

From this solution it follows that the matrix A ab c\ b in the system ( 0 ) satisfies, 

A ab c l ab h = 0 . 


(33a) 

(33b) 

(34) 


This particular result will be crucial in our later discussion. As a consequence of it, the system 
© implies another interior system on I, as no p-derivatives will arise upon evaluation on I. It 
can be solved giving, 

= -6 me 2 abcd . (35) 


a(°) 

L'abcd 


From the latter discussion that the cylinder at spatial infinity I is a characteristic of the 
system (|3(])-(31). Furthermore, because of the fact that the whole system of conformal field 
equations reduces to an interior system on I —something that does not happen with normal 
characteristics— we call it a total characteristic. 

The idea of interior systems previously discussed can be generalised by applying p times d p 
to the equations ©) and & and then evaluating on I. In this way one obtains a hierarchy of 
interior systems for the unknowns u^ = d p u\i. These quantities can be used to construct formal 
expansions of the form, 

u = Y -^-u^ p p (36) 

pi 


p —0 


for the field quantities. The resulting equations, which will be referred generically to as transport 
equations , are of the form, 


d T v (p) = Kv {p) + Q(u (0) , v (p) ) + Q(u (p) , u (0) ) 

p -1 

+ £(«< w W) jV (p-i)) +£(j)0(p-i)^ + £(p) , 


(37) 


3 =1 


which correspond to the propagation equations of the v unknowns, equations (20a)-(20g). From 
the Bianchi propagation equations (|22a|)~(|22e|) one gets, 


(V2E + A ab (c° ab )W) + A^ic^dc^ = B( 


+E 

3 =1 


(B(T^ cd y p -A - A ab y ab )^d p ^) , (38) 
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Figure 1: Spacetime close to spatial and null infinities: to the left the standard representation of 
spatial infinity as a point i°; to the right the representation where spatial infinity is envisaged as 
a cylinder. 


where C = ±. Similarly, from the Bianchi constraint equations ( 23a )-( [23q ) one obtains, 
F ab (c° ab )^d T ^ + F ab {c c ab )^dc^ p) = H{ lf b } cd , 


+ £ 

3 =1 


- F ab {cf ab )^d^-^ 


(39) 


The systems © and ( |38|) can be regarded as systems for the unknowns v^ and (f>^ if the lower 
order quantities v'-A and <f>^\ 0 < k < p — 1 are known. The two systems are decoupled from 
each other, and accordingly one would firstly solve the system (^) and then feed its solution into 
the system ( |38| ) which now could in turn be solved. 

Because of (33a) and ( 33b ), the matrix accompaigning the d T derivative in the system ( |38| ) is 
given by, 

(\/2 E + A ab (c° a J 0) ) = V / 2diag(l + r, 1,1,1,1- r) (40) 

As a consequence of this, the symbol of the system looses rank at r = ±1 —the system degenerates 
there. The points r = ±1 corresponds precisely to the sets I ± —cfr. (^) — the sets where “null 
infinity touches spatial infinity”. It is exactly this particular feature of the field equations that 
forces us to undertake a complicate and detailed analysis of the system (|3l]). From an heuristic 
point of view the degeneracy at the sets I ± can be understood as a consequence of the change 
of behaviour of the conformal boundary of the spacetime with regard to the conformal field 
equations: the cylinder at spatial infinity / is a total characteristic, while the are “only” 
normal characteristics - i.e. only subsets of the field equations reduce to interior systems on either 
or Now, standard theory of symmetric hyperbolic systems guarantees, for a given order 
p, the existence of solutions to the joint system ([37l)-(|38|) for any subset of I. However, for 


the sets I± the degeneracy implied by the equation (J40j) the usual energy estimates provide no 
information precisely at the point one is interested most—a similar phenomenon occurs with the 
original system (|30|)-(|3l|). Thus, one needs to devise non-standard methods in order to address 
existence issues —see for example the discussions in ©H- 

So far, the function k introduced in equation (jlj) has been required to be of the form k = pn', 
with k 1 analytic an such that n'(i) = 1, but otherwise it has remained unspecified [(]. Two 
choices consistent with the requirement k = pn' will be considered here. The first, n = p is 
the simplest non-trivial one. For this choice in a neighbourhood of I + is concave, while 
in a neighbourhood of I would be convex. The choice k = p has the virtue of rendering 
the simplest possible analytic expressions, both for the initial data (27a)-(27f) and the solution 
of the transport equations (|37|)-(|38|). Unfortunately, it is hard to attach to it some geometrical 

4 It is noted that the simple choice k = 1 would lead to the standard representation of spatial infinity as a point 
—see figur^ The requirement of re being of the form re = pre/ ensures that spatial infinity is blown up to a cylinder. 
See figure 111 


10 




































Figure 2: Effect of the choice of the function k on the representation of null infinity near spatial 
infinity. To the left with the choice k = p; to the right the choice n = to so that null infinity 
corresponds to the hypersurfaces r = ±1. 


significance other that its simplicity. The other choice to be considered is k = ui. This choice is 
fine as under our assumptions ui = p + 0{p). With this choice, 


0 = (1-r 2 ), 


(41) 


so that near J ± are described by the hypersurfaces r = ±1, p > 0 respectively: null infinit; 


g 


will be composed of two parallel planes, formally similar to the case of Minkowski —see 
Consequently, the system of conformal field equations @)-(|3l]) degenerate not only on / ± but 
also on Thus, the choice k = to has more geometrical and analytic relevance. As a drawback 
it renders more complicated analytic expressions. 

For future use, we note the following result on the expansion types of the diverse unknowns 
appearing in the transport equations (Et]) and (|38|) . Its proof comes from inspection O] . 


Lemma 2. The functions (c^ b — px a b)^ P \ v^\ p = 1,2,..., are of expansion type p — 2, 
p — 1, and p respectively. 


4.1 A first analysis of the transport equations 

A detailed first analysis of the structure of the solutions of the transport equations have been 
given in JIt}. Because of its relevance for the purposes of the present article, and in order to fully 
motivate our later discussion, we briefly proceed to review it. 

As seen in the previous section, the transport equations corresponding to the v ^ quantities, 
equation ( |37| ) are fully regular, thus in a first analysis one could just focus on the transport 
equations arising from the Bianchi identities, equation (|38|). The system can be written as, 


(1 + T)d T ^ + X+fiP - (p - 2)^> = Jtf, 

(42a) 

dj? + i(v + 4 p) +*-#) + # = 

(42b) 

d T ^+ l -{X + ^+X_^)=R^\ 

(42c) 

d T 4 P) + \{x+<i>P +*_#) - 4 P) = M p \ 

(42d) 

(1 - r)d T 4 P) + + (p - 2 )^ = R^. 

(42e) 


We will also make use of the transport equation arising from the Bianchi transport equations, 

5 This similarity is in some aspects deceiving, as generically when m/0 the generators of null infinity, although 
confined to the planes r = zb, are bent and may rotate a spin frame that is parallelly transported along it 
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equation ([39|). These we write as, 

Td T <t>P + X.^) -M P) = S[ p \ (43a) 

rdj p) + i(X + 4 P) - -p 4 P) = ^ p) , (43b) 

ra r 4 p) + \{X + ^ X_#) -p 4 P) = 4 P) - (43c) 

The quantities R?f \ j = 0,... ,4 and S,( p \ i = 1,2, 3, are constructed from v^ unknowns with 
0 < k < p, and from q^ quantities with 0 < k < p — 1. All these quantities are assumed to be 
known. In this first analysis their detailed structure will no be crucial. It will suffice to know 
that they are of expansion type p — 1. On the other hand, consequently with lemma |] we write, 


V 2g 

Pf= E E 

9= 12-j| fc=0 


-jj‘. : 2q.kX‘ 


2q <Z-2+j 1 


(44) 


where the coefficients a PT r. 2q .k = <ij,p-,2q,k{j) are complex valued functions. The e xpre s sion (fl4|) 
effectively introduces a decomposition into spherical harmonics for the equations (42a)-(43c). In 
this way one has only to deal with ordinary differential equations. By taking linear combinations 
of equations (42b) and (43a), ( 42c ) and ( 43b ), ( 42d ) and (13c) it can be readily seen that the 
coefficients ai^ p - 2 q ,k, d 2 ,p; 2 q,k, <i 3 ,p-, 2 q,k can be determined once the coefficients ao,p;2<2,fc and 04,p;2<2,fc 
are known. These last two coefficients can be seen to satisfy the following equations, 


(1 + T)f(T)a' 0p . 2qk + 3(r)ao,p;2<2,fe + /i(r)a 4 ,p;2g,fc = To, (45a) 

-(1 - T)/(- T )«4,p;2 g ,fe + ff(-T)a4,p;2<2,fc + ^('r)a 0 ,p;2 9 ,fc = ?4, (45b) 

where 

f(r) = 2(p+l)(p- 1) - (q- l)(q + 2)(1 - r 2 ), 

g(r) = -(p - 2)/(t) + (q- 1 )(q + 2) ^(p + 1)(1 + r) - (1 + r) 2 - ^-q(q + 1)(1 - r)(l + r) 2 ^ , 
h{r) = -?-(q - 1 )q(q + 1 ){q + 2)(1 - t) 3 , 


and To, T4, depend on the coefficients Tj'+V/'T with p < p'. They vanish, in particular, for 
q = p a s a c o nseq uence of R.‘ and S) being of expansion type p — 1. In order to solve the 
system (45a )- (|45b| ) , one can try to construct a funda men tal ma trix of the homogeneous system. 
Remarkably, the homogeneous versions of equations (45a) and (451) can be decoupled to obtain 
two Jacobi equations (second order). This class of equations have been studied extensively in the 
literature. 

For 0 < q < p — 1 one gets for the system (45a)- (|45b|) a fundamental matrix, 


X p -2q,k — 


Qp;2q,k (^”) 

( — l) 9 Qp;2<2,fc( — 7") 


where 


and 


Qp\ 2 q,k{j’') — 


Qp',2q,k (7") — 


1 — T 


1 + T 


P+2 


P-2 


(-l) q Qp-2q,k( T ) \ 

Qp;2q,k(—T) ) 

(46) 

1 

P?-f a ~ r, (r), 

(47) 

qkr 2 ' ,_2> w, 

(48) 


with P^' i5 \t) the generalised Jacobi polynomials, 

v=0 ' 


f n + a 
,n — v 


n + 0 
v 


T — 1 


(49) 
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It is noted by passing that the fundamental matrix X p - 2q ,k for given p and q does not depend on 
k. If p = q the story is other. In this case the fundamental matrix is of the form, 




p;2p,k 


CllQp-,2p,k (j~') ^12Qp-,2p,k{^~) j 

C2lQp-,2p,k{ — T) C22Qp-2p,k(—T) ) 


where, 


Qp;2p,k ) 


Qp;2p,k (l~) 




ds 

(1 + s)p- 1 (1 - s)p+ 3 ’ 


(50) 


(51) 

(52) 


and cn, C\ 2 , C21 and C22 are numerical constants. A simple induction argument shows that, 


ds 


Jo (1 + s)p- 1 (1-s)p+ 3 


A* ln(l — r) 


Ap+2 

(1 - t)p+ 2 


Ai 

(1-r) 


+ Aq 


+ 5* ln(l + t) 


Bp-2 

(1 +t)p - 2 


+ •• • + 


B 1 

(1 +r)’ 


where the A’s and the B's are some numerical coefficients. They precise values are not relevant 
for our purposes. It follows from the latter discussion that any solution for the sectors with 
p = q will contain logarithmic divergences at r = ±1 unless ao lP; 2p,fe(0) = a4, p; 2p,fc(0). Quite 
remarkably, for time symmetric initial data this quite involved and technical-looking condition 
can be reexpressed in terms of a very geometric and appealing condition on the Cotton-Bach 
tensor of the initial data at spatial infinity 0- One has the following lemma: 

Lemma 3. For time symmetric initial data the following two conditions are equivalent, 


(i) 


with p > 2 and k = 0 ,..., 2 p, 


^0,p;2p,k (0) — ^4,p;2p,/c (0) 


(53) 


(ii) 


with s = 0,1 ,... . 


D(a s b s ' ’ * ^CL\b^abcd){^) 0 ? 


(54) 


From the above discussion it should be clear that the condition (|54|) is a necessary requirement 


in order to for the solutions (iff 
to condition 


(p) 


of the transport equations (38) are smooth. Therefore we refer 
as to a regularity condition, and we have the following very suggestive result, 


Theorem 1 (Friedrich, 1997). The solutions u^ of the transport equations extend smoothly 
to the sets only if the condition 


B(a s b s ' * * B ai b 1 (i) 0, S 0, 1, . . . 


is satisfied at all orders s. It is is not satisfied for some s, the solution develops logarithmic 
singularities at I±. 


It should be emphasized that the consideration leading to this last theorem have hardly make 
use of the inhomogeneous terms of equations ( |45a| ) and (45b). The system ( 45a)-( |45b|) can be in 
matricial form as, 

Ap-,2 q (.T)y' p .2 q ,k{T) + B(T)y p . 2q ,k{T) = t p - 2q ,k , (55) 


where A p . 2q and B p . 2q are 2x2, and t p - 2q ^k is a column vector such that t p;2p ,k = 0. The the 
solution of system given in terms of the fundamental matrix is then, 


Up\2 q ,k(j) — X p . 2qt k{T)X p . 2qk (0)yp- 2qt k{Il) + 



X p-, 2 q,k( S ) A p-, 2q( S ) tp-, 2 q,k(s)ds. 


(56) 
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Assume that the regularity condition (|54|) holds. Even in this case the integrand X~^ q {r) A~ ; 2 9 (t) t p - 2 q ,k(T) 
in formula (|5(]) contains poles at r = ±1 and possibly also outside the interval [—1,1]. Therefore, 
unless some remarkable property of the conformal held equations comes into stage, the solution 
vector y P - 2 q,k will contain some logarithmic terms and consequently also the coefficients ■ 

The algebraic structure of the integral in (|b|) is too complicated to analyse directly without 
resorting to some further (still unknown) structure of the held equations. In order to shed some 
light into this direction, Friedrich & Kannar have performed an analysis of the solutions of the 
transport equations up to order p = 3 under the assumption that the initial data satishes the 
condition (^). Quite remarkably they found that no logarithms appeared in the solutions. This 
meant that somehow, the simple poles in ( |56| ) cancel out up to the order under consideration. 

At this point it is important to mention that due to the gauge conditions used in the whole 
construction, we know that the logarithmic singularities occurring in the solutions to the transport 
equations are associated with the conformal structure, and are not artifacts of a choice of gauge. 

Note that the condition equation (|54|) is conformally invariant. It allows us to identify obstructions 
to the smoothness of null infinity explicitly in terms of the initial data. Also, because of the 
hyperbolic nature of the equations suggest that if logarithmic terms are present in the solutions 
of the transport equations, then it is very likely these will propagate along . 


5 Solving the transport equations given time symmetric, 
conformally flat initial data 

On simplicity and aesthetical grounds it is natural to wonder whether the regularity condition, 
equation (|54|), is the only requirement one has to impose on the initial data in order to obtain 
solutions to the transport equations which are smooth —see for example the discussion in ]Ts| . 
Before trying to prove some statement along this lines, it is of clear interest to calculate the 
some further order in the expansions. The rationale behind it being firstly to verify whether 
the conjecture still holds, and secondly to try to find some patterns in the solutions that one 
could exploit in an eventual abstract proof. In order to simplify the calculations, we have chosen 
to restrict our attention to those time symmetric initial data sets which are conformally flat. 
These satisfy the regularity condition (|54|) in a trivial way. Thus, they represent the simplest 
(non-trivial) class of initial data sets one can look at. Their simplicity is somehow deceiving, and 
should not be regarded as a drawback on the kind of insight that can be gained through them. A 
great deal about the solutions of the Einstein field equations has been learned from the analysis 

of this class of initial data —see e.g. @,13111] • „ „ 

The already “big” transport equation systems (37)-(^8[) do not give an appropriate dimension 
of the computational difficulties one has to face if one is to take the expansions carried in to 
even higher orders. However, the calculations one has to carry are suitable to a treatment using a 
computer algebra system. In order to analyse with ease the solutions of the transport equations 
some scripts in the computer algebra system Maple V have been written. The reader is remitted 
to the appendix for further details on the implementation of the calculations on this system. 

Because of the largeness of the expressions contained inthe solutions, we have opted to provide 
a description of the qualitative features of the expansions rather than a complete list of all 
the terms calculated. In particular, attention will be focused on the solutions to the transport 
equations arising from the Bianchi identities, the functions ■ It should be emphasized that 
this does not mean that the solutions to the v unknown transport equations are not important. 
They are also crucial: a tiny mistake in the calculation of their solutions would destroy the whole 
structure of the solutions. However, as the discussion in the previous section has pointed out, the 
logarithmic terms that destroy the smoothness of null infinity appear firstly in the components 
of the Weyl spinor. 

In order to perform our expansions a number of assumptions have been made. We list them 
here for the purposes of a quick reference. 


Assumptions. It will be assumed that: 

(i) the initial data set is asymptotically Euclidean and time symmetric. 
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ii) in a neighbourhood B a (i) ofi the initial data is assumed to be conformally flat. The function 
W appearing in the conformal factor of the initial hypersurface —see equation ( \2f \ )— is a 
solution of the Laplace equation admitting in B a (i) a decomposition of the form, 


W = ^m + J2^W z p l + 0(p 9 ), 

1=1 


(57) 


where 

2 i 

W i = Y J W i ,2i,kT 2i \, (58) 

k =o 

with the coefficients Wi^i.k, i = 1,..., 7, k = 0,..., 2i complex numbers satisfying w^ 2 i,k = 
(— l) l+k Wi' 2 i, 2 i-k so that W is a real valued function. This is in consistency with the de¬ 
composition given in equation (fgjj). 

(Hi) Likewise, the components of the vector unknowns v^ and admit on I expansions in 
terms of T- \ functions consistent with lemma [|. 

(iv) The two following choices of the function n —see equation w will be considered: 

m = p, 

H2 = w. 

The result of the calculations under these assumptions are now described. 


5.1 The orders p — 0,..., 4. 

Firstly, calculations for the orders p = 1, p = 2 and p = 3 were undertaken. The results are 
in complete agreement with those given by Friedrich & Kannar when reduced to the case of 
conformally flat initial data. 

The solutions at order p = 0 can be schematically written as, 

4 0) = 0, (59a) 

<^ 0) = 0, (59b) 

= —m, (59c) 

4° } = 0, (59d) 

4 0) = 0. (59e) 


independently of the choice of n. At order p = 1 one has, 

= 0, (60a) 

<^ 1} = SX+Wfll - t) 2 , (60b) 

<j> 2 ^ = m 2 (^r 4 — 3r 2 ) + 6 W\[t 2 — 1) (60c) 

(ftp = 3X_Wi(l + r) 2 , (60d) 

4 1} = 0. (60e) 


when k = p. The expressions for (f>^\ 4>\\ <j>^ 
of <j)^ is given by, 



and for the choice k = u> are similar. That 
~1)+6W 1 (t 2 ~ 1). (60f) 
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At order p = 2 one has, 


^ 2) = fi (t)X + X+W 2 , (61a) 

<^ 2) = f2(r)mX + W 1 + f 3 (r)X + W 2 , (61b) 

<^2 2) = U( T ) m3 + h{r)mWi + / 6 (r)W 2 , (61c) 

4 2) = -^(-rK-^ - f 3 (-r)X_W 2 , (61d) 

# = fi(-r)X_X_W 2 , (61e) 

where /i(r), i = 1,... 6 are polynomials on r. Their explicit form is not relevant for our purposes. 
The polynomials are slightly different for each of the choices of k, but of the same order. Similarly, 
the components of the Weyl tensor at order p = 3 are of the form, 

4 3) = gi(r)X + X + W 3 + g 2 (r)mX + X + W 2 + g 3 (T)(X + W 7 ) 2 , (62a) 

= 94,{t)X + W 3 + g 5 {r)mX + W 2 + ge(r)WiX + Wi + g 7 (T)m 2 X + Wi, (62b) 

</4 3) = 9 s{t)W 3 + g 9 (r)mW 2 + g w (T)(Wi) 2 + gn(T) 2 m 2 Wi + g 12 (T)m A + g 13 (r)b,(62c) 
4 3) = -gi(-r)X-W 3 - g 5 (-T)mX_W 2 - g 6 (-r) WdX.Hd - g 7 (-T)m 2 X_W 7 , (62d) 
4 3) = gi (~t)X_X_W 3 + g 2 (-T)mX_X_W 2 + g 3 (-r)(X_W 1 ) 2 . (62e) 


Again, the functions gi(r), i = 1,..., 13 are polynomials, while b = 2w± i2 , 0^1,2,2 — wf 2 ± is a 
constant. 

The first new order corresponds to p = 4. Here, again, the solutions are still fully regular and 
polynomial: 


<^ 4) = /ii(t)X + ITiX + W 2 + /i 2 (t)H / iA + A + H / 2 + /r 3 (r)m(A + W'i) 2 

+ /i 4 (t)toX + A + H / 3 + h 5 (r)X + X + W 4 , (63a) 

4 4) = h e (T)iv?X+W 1 + h 7 (r)W 2 X + Wi + h 8 (r)W 7 X + W 2 + h^mW^+Wi 

+ h w (r)mX + W 3 + /in (t)X+W 4 , (63b) 

</>2 4) = h 12 (T)m 5 + h 13 (r)b + /ii^r)?™ 3 ^! + /ii5(r)WiW 2 + /ii6(t)to(Wi) 2 

+ /ii 7 (r)raW 3 + /ii8(r)W 4 , (63c) 

4 4) = -/i 6 (-t)to 3 X_1Ti - /i 7 (-t)W 2 X_1Ti - h 8 (-r)W 1 X_W 2 - /i 9 (-r)mWi.X'_Wi 

- ftio(-r)mX_W 3 - /in(—r)X_W 4 , (63d) 

= h!(-T)X_W iX_W 2 + /i 2 (-t)1TiX_X_W 2 + /i 3 (-r)m(X_TTi) 2 

+ hi(—T)mX_ X_W 3 + /i 5 (-r)X_X_W 4 . (63e) 


Again, the functions hi(r), i = 1,..., 16 are polynomials. 

5.2 The first obstructions to smoothness: orders p = 5,6. 

The calculation of the solutions to the Bianchi transport equations up to order p = 4 have shown 
that all of them are polynomial, and thus smooth at / . Consequently, the v unknowns also 
happen to be polynomial. The first modifications to this behaviour occur at the rather high order 
p = 5. Feeding the solution of the transport equations up to order p = 4 into the v transport 
equations (^) with p = 5 and solving one finds again that the components of t/ 5 ) are again 
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polynomial. However, the solution of the Bianchi transport equations are of the form, 


4 5) =C 0 m 2 G^((l-r) 7 ln(l-r) 

+ (1 + t) 3 (351 — 150r + 48 r 2 — 10r 3 + r 4 ) ln(l + r)^ + ko(r), (64a) 


4 5) = Cim 2 G^ 5) ((1 - r) 6 (2r + 5) ln(l - r) 

— (2r 3 — 15r 2 + 48r — 75)(1 + r) 4 ln(l + r)J + fci(r), (64b) 

4>^ = G 2 ?n 2 G (5) ^(1 — t) 5 (t 2 + 5t + 8) ln(l — r) 

+ (1 + t) 5 (t 2 - 5r + 8) ln(l + r)j + fc 2 (r), (64c) 

4 5) = Gi?u 2 G (5) ((1 - t) 4 (2t 3 + 15r 2 + 48r + 75) ln(l - r) 

— (1 + t) 6 (2t — 5) ln(l + r)^ + fc 3 (r), (64d) 

4 5) = G 0 m 2 G (5) ((1 - t) 3 (351 + 150r + 48r 2 + 10r 3 + r 4 ) ln(l - r) 

+ (1 + t) 7 ln(l + t)^ + fc 4 (r), (64e) 


where Go, C 4 , G 3 are non-relevant non-zero numerical factors, fci(r) i = 0 ,... ,4 are polynomials 
analogous to, for example, those in (63a)-(63e), depending on in, W\, W 2 , W 3 , W 4 , their 
derivatives and products of them. Most remarkably, 


G (5) = £Gi 5)T 4*2, (65) 

fe=0 

where, 

Gg 5) = mw 2 , 4 :,o - 2 v / 6 w 2 i2 , 0 , ( 66 a) 

Gj 5) = mw 2 ,4, 1 - 4 ^ 101 , 2 ,oWi, 2i i, (66b) 

G^ = mw2,4,2 — 4u > 2 21 — 4wi i 2 ,oWi )2 , 2 , ( 66 c) 

G^ = mw 2 ,4,3 - 4v / 3i«i, 2 ,im;i, 2 , 2 , (66d) 

G 4 5 ^ = mw 2 , 4j4 — 2 v / 6 ru 2 i2i2 . ( 66 e) 


Thus, the coefficients G^ are (besides an irrelevant numerical factor) the Newman-Penrose con¬ 
stants of the development o f th e time symmetric conformally fiat initial data — [ f 20 | . 

Plugging the solutions ( |64a|) -( |64^ ) into the transport equations for p = 6 one finds that the 
solutions of the sectors of the form T A k - develop terms containing ln 2 (l — r) and ln(l — r). 

Furthermore, the sectors T 6 k - +1 in <f>^ contain logarithms. More precisely, the solution will be 
of the form, 

^ 6) = G (5) ^lij (r) ln 2 (l - r) + Z 2 i (r) ln(l - r) + l 3j (T)\n 2 (l + r) + Lq(r)ln(l + t)) 

+ G(6) ( l 5j i T ) ln(l - r) + l 6j (t) ln(l + r)J + l 7j (r), (67) 

where j = 0,..., 7, and . (r) are polynomials. The coefficients are new obstructions to 
the smoothness of the solutions. These will be discussed a bit later. It is not hard to imagine 
that from this point onwards, terms containing ln(l — r), ln(l + r) and higher order powers of 
them will spread all around the solutions to the transport equations. Instead of analysing this 
phenomenon, we will rather focus on the smooth solutions. 
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Setting the Newman-Penrose constants to zero, the solution at order p = 6 is of the form, 


^ = D 0 m 3 G ( 6) ((2 + t)( 1 - r) 8 ln(l - r) 

+ (-254 + 233r - 128r 2 + 46r 3 - 10r 4 + r 5 )(l + r) 4 ln(l + t)) + (68a) 


^ = L»im 3 G (6) ((23 + 20r + 5r 2 )(l - r ) 7 ln(l - r) 

- (233 - 256r + 138r 2 - 40r 3 + 5r 4 )(l + r ) 5 ln(l + r)) + h{r), ( 68 b) 

<fi^ = I? 2 W 3 G < ' 6 ' ) ((64 + 69r + 30r 2 + 5r 3 )(l — t ) 6 ln(l — r) 

+ (—64 + 69r) — 30 t 2 — 5r 3 )(l + r ) 6 ln(l + r)^ + Z 2 (t), ( 68 c) 

4 6) = Dim 3 G (6) ((233 + 256r + 138r 2 + 40r 3 + 5r 4 )(l - r ) 5 ln(l - r) 

+ (23 — 20t + 5r 2 )(l + r ) 7 ln(l + r)^j + l 3 (r), ( 68 d) 

4 6) = D 0 m 3 G (6) ((254 + 233r + 128r 2 + 46r 3 + 10r 4 + r 5 )(l - r ) 4 ln(l - r) 

+ (-2 + r)(l + r ) 8 ln(l + r) j + Z 4 (r), ( 68 e) 

with lo(r),, U(t) polynomials, and 

G (6) = E G i 6)T 6 fc 3> (69) 

fc =0 

where the coefficients G^\ k = 0 ,..., 6 are given in terms of initial data quantities by, 

Gq 6) = m 2 w 3 fi,o - 12v / 10ui 3 j2 O , (70a) 

G^ = m 2 w 3 fip - 12v / Mui 2 2 0 'u;i i 2 ,i, (70b) 

G 2 6 ^ = m 2 w 3 fi ,2 ~ 24\/6wi,2,ow 2 2,1 ~ 12v / 6w 2 2 i0 w 1 , 2 , 2 , (70c) 

G^ 6) = m 2 w 3 fi 3 - 72wi i 2,o'Wi,2,iWi i2 ,2 ^ 24w 3 21 , (70d) 

G^ = m 2 W 3 , 6,4 — 24v / 6lC 2 2 j 1 1 «i, 2,2 — 12 v / 6 ’uq, 2 , 0 ^ 1 , 2,21 (70e) 

G ( 5 6) = m 2 W 3fi ,5 - 12v / 30u>1,2,1'U'i i 2,2; (70f) 

G ( 6 6) = m 2 w 3 i6i6 - 12Vl0ui 3 j2i2 . (70g) 


In analogy to the order p = 5 we will refer to these coefficients as to the order 6 Newman-Penrose 
“<constants”. One is naturally bound to ask whether the coefficients G^ are actually associated 
to some conserved quantities at null infinity in the same way that the coefficients G^ are. This 
consideration is beyond the scope of the present article, and will be analysed in detail in future 
work. It is pointed out, as a plausibility argument, that in the analysis of polyhomogeneous 
Bondi expansions carried out in for example || ^4], ^5|, the first logarithmic terms appearing 
in the expansions were associated with a conserved quantity on null infinity. As it can be seen 
from our discussion, if G^ = 0 , then the first logarithmic terms appearing in our expansions are 
precisely those associated to G^f Notwithstanding, the expansions described here are based in 
the conformal geodesics gauge, while those in ||, [2^] use the so-called Bondi gauge. Thus, one 

would have to look in detail the possible appearance of logarithmic terms in the transformation 
connecting the two gauges. 

Theorem (Main theorem, precise formulation). Necessary conditions for the development 
of initial data which is time symmetric and conformally flat in a neighbourhood B a (i ) of (spatial) 
infinity to be smooth on the set I U I + U I~ are that the Newman-Penrose constants, G\ , k = 
0,...,5, should vanish. Furthermore, the “higher order” Newman-Penrose constants, g{_ 6 \ k = 
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(5) 

0,... ,6, should also vanish. If only the coefficients G\ ; vanish then the rescaled Weyl spinor is 
at most C 2 on a neighbourhood of either I + or I~. If both and G^ vanish then generically 
the Weyl tensor will be at most of class C 3 of the latter neighbourhoods. 


From the last result one can extract directly the following (important) consequence: 

Corollary 1. The regularity condition © is not a sufficient condition for the smoothness at 
I U I + U 1“ of the development of asymptotically Euclidean, time symmetric initial data sets. 


This corollary is, thus, a negative answer to the conjecture raised in |l!|. It is nevertheless 
surprising that the obstructions to the smoothness of null infinity arise at such a high order in 
the expansions. In order to acquire a deeper understanding of why this is the case would require 
in turn an abstract understanding of the algebraic structure of the transport equations (37) and 
( |38| ) . It is not unreasonable to reckon that some group theoretical properties of the whole set up 
play a major role here. 

The theorem also suggest the following un expected conjecture: 


Conjecture. The developments of the Brill-Lindquist and Misner initial data sets possess non¬ 
smooth null infinities. 

In H H^] the Newman-Penrose constants of the Brill-Lindquist and Misner data sets [I], |l| 
have been calculated using the formula found by Friedrich & Kannar j|(|. Due to the axial 
symmetry, there is only one non-vanishing Newman-Penrose constant. Furthermore, the constant 
vanishes if and only if the data sets are actually the Schwarzschild initial data. That the devel¬ 
opment of these two initial data sets possess a non-smooth null infinity may have implications in 
the description of their late time behaviour in terms of linear perturbations on a Schwarzschild 
background. 

Some remarks regarding the theorem come also into place: 

Remark 1. The calculations for the order p = 7 are already beyond the capabilities of Maple 
V —the expressions involved are too large for the simplification routines of the computer algebra 
system, even for an Origin computer available at the Albert Einstein Institute. Nevertheless, the 
calculations of axially symmetric situations are still possible. These have been carried out for the 
orders p = 7 and p = 8 inclusive. Assuming that both G^ and G^ vanish, the solutions of the 


Bianchi transport equations have again the expected form: 

(fiP = Ejm 3 G^ (jni J (r) ln(l - r) + m2,,, (r)) + m 3j (T), (71) 

where now due to the axial symmetry there is only one order 7 Newman-Penrose constant, 

G (7) = (m 3 w 4 , 8,4 - 192^,2, 1)^8 V (72) 

If in turn G G) vanishes, 

<^ 8) = E i m 4 G (8) [n lj (r) ln(l - r) + n 2j (r)) + n 3j (r), (73) 

with, 

G (s) = (mV, 10,5 - 1920^ j2i1 )t io 5 5 . (74) 

From this evidence it is not too hard to guess the following general formula for the obstructions 
of the smoothness of null infinity in the axial situation, 

G {p) = (m p- V_3,2p-6,p-3 - 2 p_4 (p- 4)!(w 1 ,2, 1 ) p - 3 )r 2p _ 6 p_3 p _3. (75) 


The proof of such a formula is nevertheless beyond our current understanding of the transport 
equations. The significance of the latter expression will be discussed in short. 

Remark 2. In the non-axially symmetric case, it is conjectured that obstructions to the 
smoothness of null infinity (Generalised Newman-Penrose constants) are given in terms of the 
initial data by the following expressions, 

= in P Wp— 3t 2 p— 6 ,k — ] c P,k;s 0 ,Si,S 2 W l° 2 fi' w li 2 ,l W r 2 , 2 i (76) 

So + Sl+S 2 =fc 


with k = 0,1,... , p, where the coefficients c p ,fc ;S0)Sl , S2 are some numerical constants. 
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5.3 Obstructions to smoothness and the Schwarzschild initial data. 

As mentioned before, the Newman-Penrose constants of the Schwarzschild spacetime are all zero. 
Thus, in order to gain some insight into the significance of the expressions ( |75| ) and ( |76| ) it is 
convenient to see what occurs in the case of the Schwarzschild initial data. 

It is not complicate —although certainly messy— to calculate an expression for initial data for 
the Schwarzschild on the slice of time symmetry. On this slice the initial 3-metric is conformally 
flat. Thus, the required (harmonic) conformal factor can be calculated directly from the Green 
function for the three dimensional Laplace equation in spherical coordinates. The Green function 
is given by, 

00 n A — _ r 'n 

E ^^ Y n,rn^A)Y n ,m{0'A ')— T , (77) 

n —0 m=—n 

for r > r'. 

The latter expression can be lifted into the frame bundle G 0jK and written in terms of the 
functions T- k ] using formula ( |l5| ) to get, 


IT = 


00 2 n 

EE' 

n= 0 k—0 


1 rp 

-L Q 


2 n—k 


it') T, 


k n 
2 n nr 


(78) 


where {t' , r') denote the coordinates of the singularity of the Green function on the frame bundle. 
Thus, the T 2n 2n ~ k n (t') are fixed complex numbers. We write, 

Wn, 2 n,k = \mr' n n\T 2n 2n - k n (t'), (79) 

for n = 0,1,... and k = 0,..., 2 n. The latter coefficients are not independent but related to each 
other via recurrence relations, like 

(n+ 1)(2 n+ l)u>i i2 ,i x w„ )2 n,fe = \J ( k + 1)(2 n-k + l)w n+ i } 2 ( n +i),k+i 

+ r ,2 n(n+l)v / M2n _ fc)w n _i i2 (n-i),fe-i- (80) 


These can be readily obtained from similar recurrence relations holding for the spherical harmon¬ 
ics. The Schwarzschild solution has an obvious axial symmetry. If one considers an orientation 
of the coordinate system in a way that makes this symmetry of the initial data explicit —the 
singularity of the Green function is set along the z axis— one ends up with a much simplified 
expression, 

W = E(-J (ttfi.2,1 ) n T 2n n n p n , (81) 

n =0 k m / 

so that the only non-vanishing w n ^ n ,k coefficients are given by, 


^ n,2n,n 


(I) (mAi) ”" ! ■ 


1 


(82) 


One obtains precisely this expression if one requires that the hierarchy of axial Newman-Penrose 
constants, formula (f75|), vanishes at every order p. It is however important to point out that 
this is not a proof but again a plausibility argument as formula ( |75| ) has only been verified up to 
order p = 8. A much more involved calculation shows that the expressions one obtains for the 
coefficients ru 2 , 4 ,o ; • • • ,W2,aa an d w 3,6,o> ■ ■ ■ > w 3 , 6,6 by setting G^ and G^ to zero are precisely 
those one has for the Schwarzschild initial data. This involves the use of the recurrence relation 

©• 

Now, the fuction IT is a solution to AIT = 0, and thus analytic. This in turn implies that the 
function IT one deduces from requiring the coefficients ( |75| ) to vanish at all orders p is exactly 
Schwarzschildean. This evidence leads to the following, 
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Conjecture (Precise formulation). For every k > 0 there exists ap = p(k) such that the time 
evolution of an asymptotically Euclidean, time symmetric, conformal flat, conformally smooth 
initial data set admits a conformal extension to null infinity of class C k near spacelike infinity, 
if and only if the initial data set is Schwarzschildean to order p(k). 

The latter is the successor of the conjecture, in the context of time symmetric conformally 
flat, of the conjecture by Friedrich [|l9| to which reference was made in the introduction. 


6 Conclusions and extensions 


The results presented in the main theorem together with the considerations leading to the con¬ 
jecture put forward in section 5 seem to suggest that no gravitational radiation should be present 
around spatial infinity if one is to have a smooth null infinity. In other words, the notion of 
smooth null infinity seems to be incompatible with the presence of radiation around i°. Whether 
this behaviour of spacetime in the neighbourhood of spatial infinity has some implications on 
either the demeanour of the sources of the gravitational field in the infinite past or on the nature 
of incoming radiation travelling from past null infinity is a naturally arising question. 

A natural extension of the calculations here described would be to consider what happens 
with more general (i.e. not conformally flat) time symmetric initial data. If the time symmetric 
initial data is not conformally flat, then the Newman-Penrose constants are given in terms of the 
initial data by, 


q{5) 

G? 

G ( 3 5) 

py(5) 

Or4 


mw 2,4,0 
TOU>2,4,1 

mw 2 , 4,2 

■ 2,4,3 
mw 2,4,4 


2v / 6uq ; 2 i o 


254-\/6 


X 


( 2 ) 
0 ) 


4:\ /r 3wi^2,0'UJl,2,l — 


254^ 


d(2) 
ill ! 

1 


4^1,2,1 - 4u;i i 2,0«h,2 > 2 - -7 =R 


( 2 ) 

254V6“ 2 ’ 


4'\/3'U>l,2,lt/tL,2,2--pi?; 


254^6 


( 2 ) 
'3 i 


2V61U12 2- i=X 


254^6 


( 2 ) 
'4 i 


( 2 ) 

where R k ' are coefficients appearing in the expansion of the Ricci scalar r of the initial 3-metric, 
around spatial infinity. In the view of the present results one would expect to have the quantities 

~(5) 

G k appearing as obstructions to the smoothness of null infinity. Similarly, one would expect 
to obtain generalisations of the higher order Newman-Penrose constants G^ by adding 
suitable expressions containing the Ricci scalar. Ultimately, one would like to prove that there 
is an infinite hierarchy of such quantities as obstructions to the smoothness of The relation 
of these quantities with static initial data is to be analysed, at least for the first orders. It could 
well be the case that a time symmetric initial data set will yield a development with smooth null 
infinity only if it is static in a neighbourhood of spatial infinity. In relation to this, it is noted 
that the asymptotically Euclidean static initial data satisfy the regularity condition (|54|) — see 

@- 

With regard to a proof to the conjecture presented in this article, it is noted that it would 
require a much deeper understanding of the properties of the transport equations at spatial infinity 
that the one currently available. In particular one would like to be able to discuss its algebraic 
properties in a much more abstract way, and with out having to resort to “explicit expressions” 
as it was done here. As mentioned before, group theoretical properties of the set up should play 
a crucial role here. 
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Appendix: on the Computer Algebra implementation 

The calculations described in this article are hefty, however quite amenable to the use of computer 
algebra systems. In the present case the computer algebra system Maple V was used. The main 
source of difficulties in solving the transport equations implied by the conformal field equation on 
the cylinder at spatial infinity is not the actual procedure of solving the equations, but producing 
the equations that have to be solved. To understand why this is the case, we recall the quasilinear 
nature of the transport equations. Accordingly, in the construction of the equations for the diverse 
spherical sectors of the various unknowns one encounters products of the functions T-*, which 
have to be expanded using the formula ©, which in turn requires the computation of several 
Clebsch-Gordan coefficients. The number of these products heavily increases as one considers 
higher and higher expansions orders. There are Maple V scripts available in the web. We have 
made use of those written by A. B. Coates |J]. 

For the purposes of the computer algebra implementation, the diverse spinor arising have been 
considered as arrays running on the set of indices {0,1}. In this way all the elementary spinors 
introduced in section 2.1 can be readily defined. For example, the spinor x a b would be encoded 
as, 

> x[0, 0] := 0 : x[0, 1] := l/sqrt(2) : x[l, 0] := l/sqrt(2) : x[l, l] := 0 : 

The functions k j have also been considered as arrays, T[j,k,l] with empty entries, so that 
they can be manipulated formally as symbols (names, in the Maple V terminology). With this 
idea in mind, scripts that evaluate X± derivatives according to the usual rules and take complex 
conjugates of expressions containing T- have been written. Two further scripts, one implement¬ 
ing the formula (^9|) to express products of T -\'s as linear combinations of T A )’s, Texpand, and 

another that given a expression collects the terms possessing the same 7( Tcollect, are the 
core of our implementation. 

In order to compute the expansions up to a given order p one proceeds in two steps. Firstly, 
one has to compute the initial data for the conformal propagation equations. As described in the 
main text we have restricted our implementation to time symmetric, conformally flat initial data. 
The computer algebra construction of time symmetric data satisfying the regularity condition 
( fidl ) is a project on its own right that we expect to address in the near future. In the case of 
time symmetric, conformally flat data the whole free data is contained in a function W solution 
of the Laplace equation. In order to calculate the solution of the transport equations up to say, 
order p one needs to expand the function W in p up to order p. Thus, the main input of the 
script that calculates the initial data for the transport equations are the essentially the order of 
expansion p and the function kappa corresponding to the remaining conformal freedom available 
in our setting. Two choices have been used, k = p and k = u. The function W is provided as an 
object of type series in the form, 

> W := m/2 * T[0, 0, 0] + sum( , sum( , (l/(i)!) * w[i, 2 * i, k] * T[2 * i, k, i] * rho"]/, 

V = 0..2 > 1 = i)7 i' = 1..6) + 0(rho~7) : 

Two versions of the script have been written: one which deals with axially symmetric initial data, 
and the other dealing with general time symmetric, conformally flat initial data, i.e. data which 
assumes no symmetries. The computations carried in these scripts are essentially completely 
automatised (i.e. no human intervention is required), and their output is saved directly. 

The actual computation and solving of the transport equations is carried out in a second stage. 
Given an order p the scripts load first the results of the previous orders, and then proceeds to 
the most computationally expansive part of the script: the calculation of the transport equations 
at order p. As discussed in the main text, for a given order one needs first to calculate and solve 
the so-called v equations. Using our definitions of spinors as arrays, the transport equations can 
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be readily coded. For example, equation (20a) would be coded as: 


> for a in {0,1} do 

> for b in {0,1} do 

> eqnlfp, a, b] := — diff (c0[p, a, b], tau) 

— sum( , binomial(p, k) * sum( / sum( , sum( , sum( , eup[e, r] * eup[s, t] * chi[k, a, b, r, t] 

* c0[p-k,e,s]','e' = 0..1)',V = 0..1)',' s' = = 0..1)'/k' = 0..p) - f[p,a,b] : 

> od : 

> od : 


The remaining equations, (20b) to (20g), have been encoded in a similar way. The proce¬ 
dure Tcollect is then used to write all the equations in the form of linear combination of the 
T [ j ,k,l] ’s. Using, Tcollect the resulting expressions are then grouped. A further script then 
collects all the equations corresponding to a given mode and picks up the corresponding initial 
data. The systems of equations which in the generic case consists of 45 linear ordinary differential 
equations is using the dsolve procedure of Maple V with the series option set on. In this way 
Maple V calculates the solutions directly as series expansions. The order of the expansion is 
set sufficiently high, order: =60, so that the series can truncate if the solutions correspond to 
polynomials. Once the solutions have been calculated, the solutions are saved and the memory 
of the system is cleared. 

The transport equations arising from the Bianchi identities, (22a)-(22e), are likewise encoded 
and solved. As a check of correctness, the constraint equations ( |23a| )-(23c), have also been 
programmed. After the sol ution s to t he equations (22a)-( |22c| ) have been obtained, these are 
substituted in the equations (23a)-(23c) which should be satisfied automatically. As discussed in 
the main text, at order p = 5 and higher the solutions to some of the spherical sectors contain 
logarithms. These are firstly detected in the Maple V computations as series solutions that do 
not truncate. For the cases when this occurs, some supplementary scripts have been written. 
These essentially implement the analysis carried in the section 4.1. The scripts calculates the 
fundamental matrix of the system, and uses it to calculate the solution in a closed form so that 
the logarithms present in the expression are explicit. 

In order to give an impression of the progressive calculational complexity of the whole proce¬ 
dure, it is noted that the calculations in the axially symmetric case for the orders p = 1 ,... ,8 
have taken 9.6, 36, 104.2, 284.4, 894, 2917.9, 9605, and 32139 seconds respectively in a computer 
with a Dell Xeon computer (processor speed of about 2 GHz and 2Gb of RAM memory). The 
files (text only) containing the order p = 6 v-equations and b-equations are 2.2Mb and 0.7 Mb 
big respectively. 
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